Velocity and hierarchical spread of epidemic outbreaks in scale-free networks 
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We study the effect of the connectivity pattern of complex networks on the propagation dynamics 
of epidemics. The growth time scale of outbreaks is inversely proportional to the network degree 
fluctuations, signaling that epidemics spread almost instantaneously in networks with scale-free 
degree distributions. This feature is associated to an epidemic propagation that follows a precise 
hierarchical dynamics. Once the highly connected hubs are reached, the infection pervades the 
network in a progressive cascade across smaller degree classes. The present results are relevant for 
the development of adaptive containment strategies. 
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The connectivity pattern of the network of individual 
contacts has been acknowledged since long as a relevant 
factor in determining the properties of epidemic spread- 
ing phenomena This issue assumes a particular rel- 
evance in the case of networks characterized by coniplex 
topologies and very heterogeneous structures 0,0ilMlal3 
that in many cases present us with new epidemic prop- 
agation scenarios ^ ^ 0, 0, 0, 0. A striking ex- 
ample of this situation is provided by scale-free networks 
characterized by large fluctuations in the number of con- 
nections (degree) k of each vertex. This feature usually 
finds its signature in a heavy-tailed degree distribution 
with power-law behavior of the form P{k) ~ k~'^, with 
2 < 7 < 3, that implies a non-vanishing probability of 
finding vertices with very large degrees [j, Ijjl, Q • The 
latter are the "hubs" or "superspreaders" j]| responsi- 
ble for the proliferation of infected individuals whatever 
the rate of infection characterizing the epidemic, even- 
tually leading to the absence of any epidemic threshold 
below which the infection cannot initiate a major out- 
break ilfl]. This new scenario is of practical interest in 
computer virus diffusion and the spreading of diseases 
in heterogeneous populations 0, H, IE ■ 1^ also raises 
new questions on how to protect the network and find 
optimal strategies for the deployment of immunization 
resources 0, 0, ^^^^ however, studies of epi- 
demic models in complex networks have been focused on 
the stationary properties of endemic states or the final 
prevalence (number of infected individuals) of epidemics. 
The dynamical evolution of the outbreaks has been in- 
stead far less investigated and a detailed inspection of the 
temporal behavior in complex topologies is still missing. 

In this paper we intend to fill this gap by providing 
a first analysis of the time evolution of epidemic out- 
breaks in complex networks with highly heterogeneous 
connectivity patterns. We consider the time behavior of 
epidemic outbreaks and find that the growth of infected 
individuals is governed by a time scale r proportional 



to the ratio between the first and second moment of the 
network's degree distribution, r ^ (fc)/(fc^). This implies 
an instantaneous rise of the prevalence in very heteroge- 
neous networks where (fc^) — > oo in the infinite size limit. 
In particular, this result shows that scale-free networks 
with 2 < 7 < 3 exhibit, along with the lack of an intrinsic 
epidemic threshold, a virtually infinite propagation veloc- 
ity of the infection. Furthermore, we study the detailed 
propagation in time of the infection through the differ- 
ent degree classes in the population. We find a striking 
hierarchical dynamics in which the infection propagates 
via a cascade that progresses from higher to lower degree 
classes. This infection hierarchy might be used to develop 
dynamical ad-hoc strategies for network protection. 

In order to study the dynamical evolution of epidemic 
outbreaks we shall focus on the susceptible-infected (SI) 
model in which individuals can be in two discrete states, 
either susceptible or infected 0, E| ■ Each individual is 
represented by a vertex of the network and the edges are 
the connections between individuals along which the in- 
fection may spread. The total population (the size of the 
network) N is assumed to be constant and if S{t) and 
I(t) are the number of susceptible and infected individ- 
uals at time t, respectively, then N = S{t) + I{t). In 
the SI model the infection transmission is defined by the 
spreading rate A at which susceptible individuals acquire 
the infection from an infected neighbor. In this model 
infected individuals remain always infective, an approxi- 
mation that is useful to describe early epidemic stages in 
which no control measures are deployed. We shall see in 
the following that the results obtained for the SI model 
can be readily generalized to more complex schemes. 

A first analytical description of the SI model can be 
undertaken within the homogeneous mixing hypothe- 
sis 0, ^3 , consisting in a mean-field description of the 
system in which all vertices are considered as being equiv- 
alent. In this case the reaction rate equation for the av- 
erage density of infected individuals i{t) — I{t)/N (the 
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FIG. 1: a) Average density of infected individuals versus time 
in a BA network of N=10* with m = 2. The inset shows the 
exponential fit obtained in the early times (lines) and the 
numerical curves i{t) for networks with m=4,8,12,20 (from 
bottom to top), b) Measured time scale r in BA networks as 
obtained from exponential fitting versus the theoretical pre- 
diction for difi'erent values of m and N corresponding to dif- 
ferent levels of heterogeneity. 



prevalence) reads as 

di{t) 
dt 



(1) 



The above equation states that the growth rate of in- 
fected individuals is proportional to the spreading rate 
A, the density of susceptible vertices that may become 
infected, s{t) = 1 — i{t), and the number of infected in- 
dividuals in contact with any susceptible vertex. The 
homogeneous mixing hypothesis considers that this last 
term is simply the product of the number of neighbors (fc) 
and the average density i{t). Obviously, this approxima- 
tion neglects correlations among individuals and consid- 
ers that all vertices have the same number of neighbors 
(fc), i.e. it assumes a perfectly homogeneous network. At 
small times, when the density of infected vertices is very 
small, we can neglect terms of order 0{i'^) and obtain 
the leading behavior i{t) ~ ioe^^'^", where iq is the ini- 
tial density of infected individuals and th = {X{k))~^ is 
the time scale governing the growth of the infection in a 
homogeneous network. 

The above calculations recover that the outbreak's 
time scale th is inversely proportional to the epidemic 
reproduction rate; i.e. the mean number of infections 
generated by each infected individual. However, in pop- 
ulations with a heterogeneous connectivity pattern, it is 
known that the reproduction rate depends upon the con- 
tacts' fluctuations and it is natural to expect that 
also the outbreaks' time scale is analogously affected. 



Indeed, in heterogeneous networks the degree k of ver- 
tices is highly fluctuating and the average degree is not 
anymore a meaningful characterization of the network 
properties. In order to take fully into account the degree 
heterogeneity in the dynamical evolution, it is possible 
to write down the whole set of reaction rate equations 
for the average densities of infected vertices of degree k, 
ik(t) = Ik(t)/Nk, where Nk and Ik(t) are the number of 
vertices and infected vertices within each degree class k, 
respectively For the SI model the evolution equa- 

tions read as 



dikjt) 
dt 



= Xk [l^ik{t)]ek{t), 



(2) 



where the creation term is proportional to the spread- 
ing rate A, the degree k, the probabihty 1 — ik that a 
vertex with degree k is not infected, and the density Qk 
of infected neighbors of a vertex of degree k. The lat- 
ter term is thus the average probability that any given 
neighbor of a vertex of degree k is infected. In the case 
of uncorrelated networks [ijj, ©fc = © is independent of 
the degree of the emanating edge. The probability that 
each edge of a susceptible is pointing to an infected ver- 
tex of degree k' is proportional to the fraction of edges 
emanated from these vertices. By considering that at 
least one of the edge of each infected vertex is pointing 
to another infected vertex, from which the infection has 
been transmitted, one obtains 



e(t) 



j:,,{k'~i)Pik')ik>{t) 



(3) 



where (k) = J^k' k'P{k') is the proper normalization fac- 
tor dictated by the total number of edges. A reaction rate 
equation for Q{t) can be obtained from Eqs. (0) and 
In the initial epidemic stages, we neglect terms of order 
0{i^), obtaining the following set of equations 



dikjt) 
dt 

dQjt) 
dt 

These equations can 
condition ik(t = 0) 
prevalence i{t) = 



XkQ{t), 

(P) 
(fc) 



A 



1 e(t). 



(4) 
(5) 



be solved with the uniform initial 
— io yielding for the total average 
Pik)ik{t) 



i{t) = io 



where 



+ W^-W (eV--l) 



(k) 



A((fc2)_(fc))- 



(6) 



(7) 



This readily implies 
epidemic outbreak is 
Indeed, the ratio k = 



that the growth time scale of an 
related to the graph heterogeneity. 
(fc^)/(fc) is the parameter defining 
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the level of heterogeneity of the network, since the nor- 
malized degree variance can be expressed as K/{k) — 1 
and therefore high levels of fluctuations correspond to 
K ^ (fc). In networks with a Poisson degree distribution 
in which k = (k) + 1, we recover the result r ~ {X{k))^^. 
Instead, in networks with very heterogeneous connectiv- 
ity patterns k is very large and the outbreak time-scale 
T is very small, signaling a very fast diffusion of the 
infection. In particular, in scale-free networks charac- 
terized by a degree exponent 2 < 7 < 3 we have that 
AC ~ (fc^) — > 00 with the network size — > 00. There- 
fore, while T is a function of the finite size N, for large 
scale-free networks we face a virtually instantaneous rise 
of the epidemic incidence. 

It is worth stressing that these results can be eas- 
ily extended to the susceptible-infected-susceptible (SIS) 
and the susceptible-infected-removed models (SIR) jj], in 
which the Eq. ^ contains an extra term —fxikit) defin- 
ing the rate at which infected individuals of degree k re- 
cover and become again susceptible or removed from the 
population, respectively. In addition, in the SIR model 
the normalization imposes that Sk{t) = 1 — ife(i) — rk{t), 
where rk{t) is the density of removed individuals of de- 
gree k. The inclusion of the decaying term —fiik, does 
not change the picture obtained in the SI model. By us- 
ing the same approximations, the time scale is found to 
behave as r ~ (fc)/(A(fc^) — {n + l)(fc))- In the case of di- 
verging fluctuations the time-scale behavior is therefore 
still dominated by {k'^) and r is always positive whatever 
the spreading rate A. This allows to recover the absence 
of an epidemic threshold, i.e. the lack of a decreasing 
prevalence region in the parameter space. 

In order to check these analytical results we have per- 
formed numerical simulations of the SI model on net- 
works generated with the Barabasi-Albert (BA) algo- 
rithm 2^. These networks have a scale- free degree dis- 
tribution P(fc) k~^ and we use different network sizes 
TV and minimum degree values m in order to change 
the level of heterogeneity, that in this case is given by 
K m\nN. Simulations use an agent-based modeling 
strategy in which at each time step the SI dynamics is ap- 
plied to each vertex by considering the actual state of the 
vertex and its neighbors. It is then possible to measure 
the evolution of the average number of infected individ- 
uals and other quantities, by averaging over 10^ realiza- 
tions of the dynamics. In addition, given the stochastic 
nature of the model, different initial conditions and net- 
works realizations can be used to obtain averaged quan- 
tities. In Fig. n we report the early time behavior of 
outbreaks in networks with different heterogeneity lev- 
els and the behavior of the measured r with respect to 
(fc)/A((fc^) — (fc)). The numerical results recover the an- 
alytical prediction with great accuracy. Indeed, the BA 
network is a good example of heterogeneous network in 
which the approximations used in the calculations are 
satisfied [l9| . In networks with correlations we expect to 




FIG. 2: a) Time behavior of the average degree of the newly 
infected nodes for SI outbreaks in BA networks of size A'' — 



10" 



Time is rescaled by r. Reference lines are drawn at 



the asymptotic values {k^) / {k) for t <^ r and m for t ^ t. 
The two curves are for m — 4 (bottom) and m — 14 (top), 
b) Inverse participation ratio Y2 versus time for BA network 
of size N — 10'' with minimum degree m = 4,6, 8, 10, 12, 14 
and 20, from top to bottom. Time is rescaled with r. The 
reference line indicates the minimum of Y2 around t/r ~ 6.5. 



find different quantitative results but a qualitatively sim- 
ilar framework as it happens in the case of the epidemic 
threshold evaluation 19]. 

The previous results show that the heterogeneity of 
scale-free connectivity patterns favors epidemic spread- 
ing not only by suppressing the epidemic threshold, but 
also by accelerating the virus propagation in the pop- 
ulation. The velocity of the spreading leaves us with 
very short response times in the deployment of control 
measures and a detailed knowledge of the way epidemics 
propagate through the network could be very valuable 
in the definition of adaptive strategies. Indeed, the epi- 
demic diffusion is far from homogeneous. The simple for- 
mal integration of Eq. |(2Jl written for Sk = I — ik yields 
Sk{t) = sOe-^'=*(*) where <i>{t) = J^dt'e{t'). This resuh 
is valid for any value of the degree k and the function 
^(t) is positive and monotonically increasing. This last 
fact implies that Sk is decreasing monotonically towards 
zero when time grows. For any two values k > k' , and 
whatever the initial conditions s° and s°, are, there ex- 
ists a time t* after which suit) < Sk'it). This means that 
vertices belonging to higher degree classes are generally 
infected more quickly. A more precise characterization 
of the epidemic diffusion through the network can be 
achieved by studying some convenient quantities in nu- 
merical spreading experiments in BA networks. First, we 
measure the average degree of the newly infected nodes 
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at time t, defined as 



hnfit) 



Efefc[4W-/fc(t-l)] 
I{t)-I{t-1) ■ 



(8) 



In Fig. Oi) we plot this quantity for BA networks as a 
function of the rescaled time t/r. The curves show an 
initial plateau that can be easily understood by consid- 
ering that at very low density of infected individuals i, 
each vertex will infect a fraction of its neighbors without 
correlations with the spreading from other vertices. In 
this case each edge points to a vertex with degree k with 
probability kP{k)/{k) and the average degree of newly 
infected vertices is given by fcinf(i) = (fc^)/(fc). After 
this initial regime, fcinf (t) decreases smoothly when time 
increases. The dynamical spreading process is therefore 
clear; after the hubs are very quickly infected the spread 
is going always towards smaller values of k. This is con- 
firmed by the large time regime that settles in a plateau 
kinf{t) = m; i.e. the vertices with the lowest degree are 
the last to be infected. 

Further information on the infection propa gati on is 
provided by the inverse participation ratio Y2 {t) [21| . We 
first define the weight of infected individuals in each class 
degree k by Wk{t) = Ik{t)/I{t). The quantity Y2 is then 
defined as 



(9) 



If Y2 ~ l/fcmax, infected vertices are homogeneously dis- 
tributed among all degree classes. In contrast, if Y2 is 
not small (of order unity) then the infection is localized 
on some specific degree classes that dominate the sum 
of Eq. @. In Fig. 133) we report the behavior of Y2 ver- 
sus time for BA networks with different minimum degree. 
The function Y2 has a maximum at the early time stage, 
indicating that the infection is localized on the large k 
classes, as we infer from the plot of see Fig.|2t)- 

Afterwards Y2 decreases, with the infection progressively 
invading lower degree classes, and providing a more ho- 
mogeneous diffusion of infected vertices in the various k 
classes. Finally, the last stage of the process corresponds 
to the capillary invasion of the lowest degree classes which 
have a larger number of vertices and thus provide a larger 
weight. In the very large time limit, when the whole net- 
work is infected, Y2{t = 00) = J2k^i^)^- Noticeably, 
curves for different levels of heterogeneity have the same 
time profile in the rescaled variable t/r. This implies 
that, despite the various approximations used in the cal- 
culations, the whole spreading process is dominated by 
the time-scale defined in the early exponential regime of 
the outbreak. 

The presented results provide a clear picture of the in- 
fection propagation in heterogeneous networks. First the 
infection takes control of the large degree vertices in the 



network. Then it rapidly invades the network via a cas- 
cade through progressively smaller degree classes. The 
dynamical structure of the spreading is therefore char- 
acterized by a hierarchical cascade from hubs to inter- 
mediate k and finally to small k classes. This result, 
along with the very fast growth rate of epidemic out- 
breaks, could be of practical importance in the set-up 
of dynamic control strategies in populations with het- 
erogeneous connectivity patterns. In particular, targeted 
immunization strategies that evolve with time might be 
particularly effective in epidemics control. 
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